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An analysis of the properties of half-doped manganites is presented. We build up the phase 
diagram of the system combining a realistic calculation of the electronic properties and a mean 
field treatment of the temperature effects. The electronic structure of the manganites are described 
with a double exchange model with cooperative Jahn- Teller phonons and antiferromagnetic coupling 
between the Mn core spins. At zero temperature a variety of electronic phases as ferromagnetic 
(FM) charge ordered (CO) orbital ordered (00), CE-CO-OO and FM metallic, are obtained. By 
raising the temperature the CE-CO-OO phase becomes paramagnetic (PM), but depending on 
the electron-phonon coupling and the exchange coupling the transition can be direct or trough 
intermediate states: a FM disorder metallic, a PM-CO-00 or a FM-CO-OO. We also discus the 
nature of the high temperature PM phase in the regime of finite electron phonon coupling. In this 
regime half of the oxygen octahedra surrounding the Mn ions are distorted. In the weak coupling 
regime the octahedra are slightly deformed and only trap a small amount of electronic charge, 
rendering the system metallic consequentially. However in the strong coupling regime the octahedra 
are strongly distorted, the charge is fully localized in polarons and the system is insulator. 

PACS numbers: 75.47.Gk,75.10.-b. 75.30Kz, 75.50.Ee. 



I. INTRODUCTION 

Oxides of composition (Ri_a;Aa;)Mn03 where R de- 
notes rare earth ions and A is a divalent alkaline ion, 
are called generically manganites. In these compounds 
X coincides with the concentration of holes moving in 
the eg orbital band of the Mn ions that ideally form 
a cubic structure. The electronic and magnetic prop- 
erties of manganites are determined by the competition 
between at least four independent energy scales; the anti- 
ferromagnetic interaction between the Mn spins, the elec- 
tron phonon coupling, the electronic repulsion and the 
kinetic energy of the carriers. In manganites the en- 
ergy magnitude of these effects is the same and very 
different states can have very similar energies. Conse- 
quently slightly varying parameters as carrier concentra- 
tion, strain, disorder, temperature etc, different phases 
can be experimentally observed. 

In this work we calculate the phase diagram of man- 
ganites at half doping, x=l/2. We study the different 
phases of the system as function of temperature (T), 
electron-phonon coupling (A) and antiferromagnetic su- 
perexchange coupling between the classical Mn t2g core 
spins (Jaf)- In particular we are interested in analyzing 
the different phase transitions that undergoes the para- 
magnetic phase when the temperature is lowered (Ij, 

Previous temperature phase diagrams have been ob- 
tained for the one orbital double exchange model'^,'^,'^. 
Although this model clarify some physics of the mangan- 
ites at low hole doping, the one orbital model is insuf- 
ficient for describing the complexity of half doped man- 
ganites. On the other hand Monte Carlo simulation of 
a two orbital double exchange model are very expen- 
sive computationally and only very small systems can 



be studiedly]. In this paper the phase diagram is ob- 
tained combining the electronic properties obtained from 
a realistic microscopic two-dimensional Hamiltonian and 
a mean field treatment of the temperature fluctuations. 
With this method is possible to analyze system with up 
to 26x26 Mn ions, and the results are rather free of finite 
size effects. It has been shown extensibelyQ that, due 
to the directionality of the electronic active Cg Mn or- 
bitals, the phases appearing in two-dimensional calcula- 
tions are qualitatively similar to those appearing in three- 
dimensional models, and the zero temperature phase di- 
agrams obtained in two and three dimensions are topo- 
logically equivalents. 

In Figures I1I2I and |21 we show the main results of 
this work. Fig^ shows the phase diagram T — A of 
the half-doped two-dimensional, two orbital double ex- 
change model with cooperative Jahn- Teller phonons for 
JAF=0-06t In FigEl and we plot the T - Jaf phase 
diagrams for \~1.6t and X~2t respectively. Here t is the 
hopping between first neighbors dx^-y^ orbitals in the x- 
y plane, and in the rest of the paper it is taken as the 
unit of energy. 

The magnetic phases appearing in the diagram are: 
ferromagnetic (FM), antiferromagnetic (AF), paramag- 
netic (PM) and CE (ferromagnetic zigzag chains that are 
coupled antiferromagnetically0). In some phases the 
electric charge is chessboard-like ordered (CO phases). 
The charge modulation is due to the cooperative Jahn- 
Teller coupling. The electronic charge trapped in a AIn 
ion depends on the distortion of the oxygen octahedron 
surrounding the ion, that at the same time depends on A. 
The octahedra distortions of first neighbors Mn ions are 
coupled and due to the elastic energy a distorted octa- 
hedron is encircled by undistorted octahedra. The order 
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FIG. 1: Phase diagram, T-X, for the two dimensional DE 
two orbital model with cooperative Jahn- Teller phonons. The 
AF coupling is JAF=0.Q6t. The abbreviations explaining the 
phases are defined in the text. First order transitions are 
indicated by continuous lines and second order transitions by 
dashed lines. The big dote line indicates a metal-insulator 
transition. The PM-M is also CD, and the FM-CO is also 
OO. Note that for lAt < A <1.5t and low T there is a smaU 
region, almost imperceptible in the figure, where the system 
is FM-CO-OO. 
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FIG. 2: Phase diagram, T-Jaf, for the two dimensional DE 
two orbital model with cooperative Jahn- Teller phonons. The 
Jahn- Teller coupling is A=1.6t. The abbreviations explaining 
the phases are defined in the text. First order transition are 
indicated by continuous lines and second order transitions by 
dashed lines. 



FIG. 3: Same than Fig. El but with \=2t. 



in the octahedra distortions produces the electron charge 
modulation. For large values of the coupling, A > 2.3t, 
the charge modulation is maximum and the CO phase is 
an ordered distribution of Mv?^ and Mn^^ . However for 
weaker values of A the modulation is much smaller. The 
cooperative nature of the Jahn- Teller coupling implies 
that the CO is concomitant to the existence of orbital 
ordering (00). Charge disordered (CD) names metallic 
phases (M) where the electric charge moves freely be- 
tween the Mn ions and is equally distributed between 
them. For strong Jahn- Teller coupling, A > 1.9i, by rais- 
ing T the CO phase melts into an insulator polaron-like 
state where the carriers are practically fully localized in 
M n ions surrounded by a strong distorted oxygen octa- 
hedron. For weak values of A the CO phase melts into 
a state formed of weak deformed oxygen octahedra, that 
trap only a small amount of the electronic charge. This 
phase is metallic and can be considered as a crossover 
between the normal Fermi liquid phase and the polaron- 
like state|^0. In the phase diagrams, Figures [TO and 
O the solid and the dashed lines represent first and sec- 
ond order phase transitions respectively. In Fig ^ the 
large dot line separating the disordered octahedra metal- 
lic phase from the polaron-like state, is an estimation of 
the value of A where the system undergoes a metal in- 
sulator transition. Some of the high temperature phases 
appearing in Figures 1 1121 and have been not reported 
in Monte Carlo simulations. We believe that due to size 
effects the disordered octahedra phases could be difficult 
to identified in small systems. 

Due to computational limitations we have only con- 
sidered spatially uniform phases in Figures I1I2I and O 
. Therefore we can not rule out the existence of more 
sophisticated spatiall y commensuarte /inconmensurate 
modulated phaseslH 111 El 111 III El • In the calcu- 
lation we have not considered coupling between different 
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order parameters. This could affect or phase diagram 
near crossing between different lines of phase transitions. 
Also near the crossing points the thermal fluctuations 
that are neglected in the mean field approximation could 
be important. 

Some characteristics of the phase diagrams explain sev- 
eral experimental results on manganites and deserves to 
be discussed explicitly: 

i) There is a region on the parameter space, small Jaf 
large A, where the zero temperature ground state is FM 
and presents charge and orbital ordering, FM-CO-00. 
This phase results from the cooperative JT phonon and 
does not requires a strong Jaf couphng to exists [^ll^. 
A phase with the same symmetry has been observed re- 
cently by Loudon et aLjl8|. 

ii) I n ag reement with several experimental results 1 191 1201 
111 111 in El and theoretical works 0, IM El IH Ell 

there is a wide region of the phase diagram (intermediate 
values of A and Jaf)i where the low temperature ground 
state is the insulator CE-CO-00 phase. 

iii) At high temperature the system is always paramag- 
netic, but depending on the value of the electron phonon 
couphng the system behaves as a metal(low A) a disor- 
der metal (moderates values of A) or an insulator (strong 
coupling). These regimes have been observed experi- 
mentally. In the case of Lao.5Sro.5Mn03(29j the high 
temperature paramagnetic phase is metallic, whereas 
in perovskites with strong electron-phonon coupling, as 
Lao.sCao.sMnOsP] the paramagnetic phase is an insula- 
tor. 

iv) By raising the temperature the CE-CO-00 insulating 
phase becomes a PM phase, but depending on the values 
of A and Jaf, in the heating process the transition is 
trough a FM disorder metallic phase or through a PM- 
CO-00 phase. These phase transition paths describe 
appropriately the experimental phenomenology. Upon 
coolingjLao.sCao.sMnOa first becomes FM-CO and then 
CEjll |30|, however in the case of Pro, 5 Cap. 5 MnO 3, the 
intermediate phase is a PM-CO phaselalia 

This article is organized as follows. In Sec. II we de- 
scribe the model and the approximations used to find 
the energies and properties of the different phases. In 
Sec. Ill, we describe the phases studied and outline the 
method for obtaining the critical temperatures and the 
transitions between the phases. In Sec. IV we summarize 
the results. 



II. MODEL. 

We assume that the manganites crystallize in an ideal 
perovskites structure where the Mn ions form a cubic 
lattice. The crystal field splits the Mn d levels into an 
occupied strongly localised t2g triplet and a doublet of eg 
symmetry. The Coulomb interaction between electrons 
prevents double occupancy and aligns the spins of the d 
orbitals. At x ^1 there is a finite density of electrons in 
the system that hop between the empty Cg Mn states. 



The Hund's coupling between the spins of the carriers 
and each core spin is much larger than any other energy 
in the system, and each electron spin is forced to align lo- 
cally with the core spin texture. Then the carriers can be 
treated as spinless particles and the hopping amplitude 
between two Mn ions is modulated by the spin reduction 
factor, 

fi2 = COS Y COS y -f e*^*'! sm y sm y (1) 

where (j>i} are the Eulcr angles of the, assumed clas- 
sical, Mn core spins {S^} . This is the so called DE 

model m mm. 

For obtaining the states of the system we study a 
two dimensional DE model coupled to Jahn- Teller (JT) 
phonons. We also include the AF coupling between the 
Mn core spins Jaf- 

i,j,a,a' 

+ Jaf J2 S.S, + i^(/3Q?,-f gi. + QiO 

<i.j> i 

+ ^'^{QliPi+ Q2iTxi+ QziTzi) , (2) 

i 

here Cf^ creates an electron in the Mn ions located at 
site i in the Cg orbital a (a=l,2 with \—\x^ — > and 
2= 1 32;^ — >). The hopping amplitude is finite for next 
neighbors Mn and depends both on the type of orbital 
involved and on the direction u between sites i and j 
{tlf = ±V?>tlf = ±Vitlf = itf:^'^ =t)%t\s taken 
as the energy unit. The forth term couples the eg elec- 
trons with the three active MnOg octahedra distortions: 
the breathing mode Qu, and the JT modes Q2i and 
Q-ii that have symmetry x^-y^ and Zz^-r^ respectively. 
Q\i couples with the charge at site i, pi = Cf^Ci^a 
whereas Q2i and Qsi couple with the x and z orbital pseu- 
dospin, Txi = C^Ci2 + C^^Ca and t^_, = C^^Cn - C%Ci2, 
respectively. The third term is the elastic energy of the 
octahedra distortions, being /3 >2 the spring constant 
ratio for breathing and JT-modes0. In the perovskite 
structures the oxygens are shared by neighbouring MnOg 
octahedra and the Q's distortions are not independent, 
cooperative effects being very important 0. In order to 
consider these collective effects we consider the position 
of the oxygen atoms as the independent variables of the 
JT distortions. 

For a given value of the parameters A and Jaf, and 
a texture of core spins {S^}, we solve self-consistently 
the mean field version of Hamiltonian (01 and obtain the 
energy, the local charges {pi}, the orbital pseudospin or- 
der {TxijTzi} and the oxygen octahedra distortions Qa^i- 
These quantities are better described by their Fourier 
transforms, that are represented by the same symbol with 
a hat: /5(G), (3i(G), ... We have checked that with our 
model and the appropriated parameters we recover some 
well known states as the FM and the CE-COHHEil. 
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FIG. 4: Density of states for different ordered phases (see 
text). Tfie dasfied lines indicate the position of the Fermi 
energy. The density of states of the FM-CO-OO phase is 
discrete as it is obtained by solving numerically a system with 
24x24 Mn ions. 



Therefore we are in a position for studying other phases 
that can appear at half doping. 



III. ELECTRIC AND MAGNETIC PHASES 
AND CRITICAL TEMPERATURES. 

A. Ferromagnetic Metallic (FM-M) phase. 

In this phase the Mn spins point, on average, in a 
particular direction, and there is a finite relative mag- 
netization < m >. In the spirit of the virtual crystal 
approximation, we consider a unique value for the spin 
reduction factor that corresponds to its expectation 
value. 



' 1 + COS 9i 



1+ < COS 9ij > 



where < > means thermal average, 0y represents the an- 
gle formed by the spins located at sites i and j. In Eq.© 
we have approached the expectation value of the square 
root by the square root of the average value's?]. This 
phase occurs for small values of Jaf and A for which the 
octahedra are not distorted and the system is metallic. 
In Fig^a) we plot the density of states (DOS) of the 
FM metallic phase for < m >=1. For other values of the 

magnetization the bandwidth is 8t^ i+^n>f._ r^^ie inter- 
nal energy is the sum of an electronic energy and the AF 
superexchange energy, 



^ — £fm 



1+ < m >2 



2Jaf <m>^ 



(4) 



Here Sf m is the electronic energy per Mn ion in the full 
polarized (< m >=1) ferromagnetic metallic phase. The 
paramagnetic phase corresponds to a random orientation 
of the AIn spins, < m >=0. Both the ferromagnetic 
and paramagnetic metallic phases do not present charge 
ordering, however in order to minimize the kinetic en- 
ergy ,the x^-y^ orbital is more occupied than the Sz^-r^ 
orbitalj^ and these phases show orbital ordering. The 
occurrence of orbital order in a metallic phase only hap- 
pen in two dimensions and therefore, here, we do not 
consider it as an relevant effect. 

In order to describe thermal effects we have to com- 
pute the free energy. In the manganites, the typical 
critical temperatures are lower than 300i^ and the ther- 
mal energy is much smaller than the electronic energies, 
t ^--^0.26^, and neglect thermal effects on the carriers 
Fermi distribution is a good approximation. Therefore 
we only need to consider the entropy of the classical Mn 
spins. We use the molecular field approximation, neglect- 
ing all correlations between different spins and assuming 
for each of them a statistical distribution corresponding 
to an effective magnetic field /i|3^. The free energy of a 
system of classical spins in presence of h is 



and the magnetization of this system is 

dT 1 1 



<m>= — 



dh tanh(/3/i) f3h ' 



(5) 



(6) 



This equation gives the relation between the effective 
magnetic field and the magnetization. The entropy of 
the spin system is given by 



-TS = T- <m> h 



log 2 



. sinh(/3/i) \ 
/3h J 



m,j3h 



(7) 

that can be written as a function of < m > solving self- 
consistently Eq^l and EqEl In the limit of small magne- 
tization the entropy takes the form. 



(3) 



S{<m>) 



log 2 
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FIG. 5: (Color online) Schematic view of CO-00 phases. 
The horizontal and vertical lines on the lattice sites represent 
the positive and negative oxygen octahedra distortions. The 
size of the dots is proportional to the charge located on the 
Mn ions. In the case of A=1.6f the octahedra distortion is 
Q2~H-25t and the ions have charge 0.71 and 0.29. In the 
FM phase all ions have the same spin orientation, whereas 
the CE phase consists of FM zigzag chains which are coupled 
AF. 

Using this expression for the Mn spins entropy the total 
free energy of the system near the PM to FM transition 
takes the form, 

F{<m>) = E^^^-^^ -TS{<m>) 

« F{0)+<m>^i^T+'^ + 2JAF) 

+ <m>^dT-'-^) + ... (9) 
20 8\/2 

and the Curie temperature is 

The constant muhiplying the <m>^ term is positive and 
therefore the transition is a second order phase transition. 



charge and orbital order and a finite relative magneti- 
zation, < m >y^ 0. This is a surprising state as the 
charge order (CO) is normally associated with AF order. 
In manganites the CO is created by the electron phonon 
coupling and for realistic values of A, the charge dispro- 
portionation between the two different types of Mn ions 
is small. Then the system can reduce its energy by order- 
ing ferromagnetically the Mn core spins and minimizing 
the kinetic energy. 

The charge ordering is characterized by a single Fourier 
component, (tTjTt), of the charge density, /5(7r,7r) ^ 0. 
The orbital ordering is described by a uniform z-orbital 
pseudospin, Tz^i and two finite Fourier components of 
the a;-component of the pseudospin, r2;(7r/2, 7r/2) = 
Ta;(— 7r/2, — 7r/2) ^ 0. Due to the electron phonon 
coupling, the amplitude of the charge modulation is 
related with the distortion of the oxygen octahedra, 
Q2(7r/2,7r/2) = Q2{-Tr/2,-n/2) ^ 0, that depends on 
the value of A. The sites where the electron charge is 
larger than the average, a:=l/2, the Q2 mode has a finite 
value, whereas the octahedra surrounding Mn ions with 
less charge than the average are undistorted. The sign 
and amplitude of the distortion are modulated in order 
to minimize the elastic energy of the cooperative Jahn- 
Teller distortion. The CO-00 phase can be described as 
a distribution of deformed octahedra. The position of the 
distorted octahedra determines the electron charge order- 
ing, and the orientation of the octahedra distortions fix 
the orbital ordering in the system. CO and 00 are a 
consequence of the spatial modulation of the octahedra 
deformations. In Fig. \E\we plot the arrangement of the 
octahedra distortion and the electronic charge distribu- 
tion for a unit cell with 8x8 Mn ions for X=1.6t. As the 
Mn ions have two orbital per site, the (7r/2,7r/2) mod- 
ulation of the octahedra distortion opens a gap at the 
Fermi energy for a;=l/2 and therefore the system is an 
insulator. The magnitude of the gap increases with the 
value of A. In Fig^Jb) we plot a typical density of this 
phase for A=1.6<. 

A decrease of the relative magnetization < m > pro- 
duces a reduction in the spin reduction factor fij and in 
the kinetic energy, making the relative electron-phonon 
interaction stronger. The internal energy per Mn ion of 
this phase can be written as 

-eFA/~co(A,<m>) + 2JAF <m>2 . (11) 

where the electronic energy depends in a complicated way 
on A and < m > and has to be obtained numerically by 
solving Eq. Q . 

The Curie temperature of the FM-CO-00 phase can 
be calculated similarly to the FM-M case, and we obtain 



B. Ferromagnetic, Charge Ordered, Orbital 
Ordered phase (FM-CO-OO). 

For small values of Jj^p and strong enough electron 
phonon coupling, A >1.4t, the ground state present 



Tc = -- 



2 a£;^M-co 



3 9 <m>2 



(12) 



<m>=0 



In the case of the FM-CO-00 phase the derivative has 
to be obtained numerically. From higher derivatives of 



6 



the internal energy with respect to the magnetization we 
obtain that the transition is second order. 



C. CE, Charge Order, Orbital Order phase 
(CE-CO-OO). 

MagneticaUy the CE phase consists of ferromagnetic 
zigzag chains coupled ferromagnetically. The horizontal 
and vertical steps of the chain contain three Mn ions. For 
finite values of the electron phonon coupling the charge, 
orbitals and octahedra distortions are ordered similarly 
to in the FM-CO-00 phase. In the case of A=0, the 
CE phase does not present charge order, however due 
to the dependence of the electron hopping in the spatial 
direction, the system spontaneously creates an orbital or- 
der similar to that occurring in the FM-CO-00 phase, 
T,(0,0) ^ and T,(f,f) = Tx(-f,-f) 7^ 0. Indepen- 
dently of the existence of charge modulation the (f , f ) 
orbital order opens a gap at the Fermi energy an the half 
doped system is an insulator. In the A=0 case the gap is 
just due to the orbital order and the system is classified 
as a band insulator ;35] . In Fig. 01c)and(d), we plot the 
density of states in the CE phase for X=1.6t and A=0 
respectively. There is a gap at the Fermi energy that for 
A=0 is of the order of the hopping parameter, t, and that 
increases with the electron phonon coupling. 

For each value of A there is a minimal value of Jaf 
for the occurrence the CE-OO-CO phase. At zero tem- 
perature and for A=0 the minimal AF coupling is Jaf — 
O.lSt, whereas for A=1.6i and X=2t the minimal AF cou- 
pling are 0.05t and 0.023t respectively. Figures |21 and 01 

The magnetization of the CE phase is described by 
the relative amount of saturation in each zigzag chain 
< ms >■ In the virtual crystal approximation the fluctu- 
ations are neglected and the hopping is modulated by the 
spin reduction factor that is different along the zigzag FM 
chain, f^^^ than between the AF coupled chains. f^^[33|. 



rFM 



f 



AF 



1+ <r7l5>2 



1- <ms>2 



(13) 



The internal energy of this phase depends on A, and < 
ms > and can be written as. 



E 



CE-CO 



ecE-co{\<ms>) 



(14) 



As each Mn spin core is surrounded by two Mn spins 
coupled FM and other two coupled AF, the superex- 
change energy is zero. 

In order to compute the Neel temperature of the CE 
phases, we introduce an effective field for each spin sub- 
lattice. Taking into account that both, the magnetiza- 
tion and the effective magnetic field, have different sign 
in each sublattice, we obtain the same expression for the 
entropy than in the FM phase, but just changing < m > 



by < mg > 33] . With this the free energy takes the form, 
F(<ms>,A) = E^^-^o -TS{<ms>) 



F{0)+ <ms>^ (nT + a 



(15) 



with 



^CE-oo oicte + a <ms>^ +b<ms>^ +.... (16) 

Due to the symmetry of the CE phase the coefficient a is 
zero, and the Neel temperature gets the form 



t; 



CE 
N 



20 

I 

9 



(17) 



Numerically, we obtain 6 < and therefore there is a 
finite temperature first order phase transition. As it is 
shown in Fig. I1I3I only for large values of A and Jaf the 
CE-00-00 phase transforms directly into the PM-CO- 
00. At smaller values of A and Jaf the transition to the 
PM phase is through a FM-CO-00 (FigUl or through a 
FM metallic disordered octahedra phase (Fig. [2J). 



D. Antiferromagnetic CO-OO phase 

For large values of Jaf, the kinetic energy is zero 
and the electron charge is fully localized forming a 
chessboard-like pattern of Mn^^ and Mn'*+. For finite 
values of A the octahedra surrounding the Mn'^+ ions are 
distorted cooperatively and the system presents 00 and 
has the same symmetry than the others CO-OO phases. 
The minimal value of Jaf for the occurrence of this phase 
depends on A, but in general is rather large, and this 
phase is unlikely to occur in half doped manganites, and 
it in this work is not considered. 



E. Disordered octahedra phases: Polaron phase 
and metalHc crossover phase. 

At high temperature the CO phases melt. The nature 
of the disorder high temperature phase depends on the 
strength of the electron phonon interaction. The melt- 
ing of a charge density wave is a complicated problem, 
and only recently it has been possible to apply an unique 
treatment for studying the weak and the strong coupling 
regimes 9, 10]. Part of the difficulty in studying the melt- 
ing of a charge density wave is to define the order param- 
eter that vanishes in the disordered phase. In the man- 
ganites this problem can be avoided by considering the 
distribution of the distorted oxygen octahedra instead of 
the modulation of the charge density. In the CO-OO 
phases there are two order parameters, ^, describing the 
order-disorder transition in the the spatial distribution 
of the distorted oxygen octahedra and rj, parameterizing 
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FIG. 6: (Color online) Schematic view of the disordered oc- 
tahedra phase for moderate values of A. The horizontal and 
vertical lines on the lattice sites represent the positive and 
negative oxygen octahedra distortions. The size of the dots 
is proportional to the charge located on the Mn ions. For 
\=1.6t the octahedra distortion is Q2=±0.95t, and the charge 
located on the Mn ions is in the range 0.25-0.75. 



their orientational order. The octahedra distorted phase 
can exists on a FM or PM Mn spins background. How- 
ever, we have not found any region in the parameter space 
where the stable phase contains a disorder distribution of 
octahedra on a magnetic CE background. 

The free energy of the CO-00 phases has the form, 

= F(0,0)-T (sm{0 + ^^A/(r/)) r?) 



being Sm the entropy of mixing [ 



(18) 



Sm{x) ^\n2-^{l + x)\n{l + x) - ^(1 - x) ln(l - x) , 

(19) 

and E'-''^ the electronic internal energy mea- 

sured with respect to the full disordered case. Near ^=0 
and 77=0, the internal energy can be expanded as 



E 



CO-00 



(20) 



we have evaluated numerically the lower coefficients in 
this expansion and we have found a strong coupling be- 
tween the two order parameters that indicates the first 
order nature of the order disorder transition. Assuming 
a transition between the fully disordered phase and the 
fully ordered phase the melting critical temperature of 



the CO-00 phases is 
2 



Tm — 



31n2 



E- 



CO-00 



(1,1) -E 



CO-00 



(0,0)) (21) 



In order to compute the energy of the disordered octa- 
hedra phase we locate randomly the distorted octahedra 
(^=0). The sign of the distortion is also random (?7=0). 
We assume that all the octahedra have the same elonga- 
tion amplitude, Q2, and we find its value by minimizing 
the energy. For a given value of the electron phonon cou- 
pling the value of the deformation Q2 is smaller in the dis- 
order phase than in the CO-00 phases. We have checked 
in several cases that the energy is minimized when the 
number of distorted octahedra and the number of elec- 
trons is the same. 

By changing the value of A, the octahedra disordered 
phase has different behaviors. There is a critical value 
of A where the distortion Q2 become finite. As can be 
observed in Fig^ this value, A ^1.3t, is slightly smaller 
than the critical A for the zero temperature FM to CO- 
00 transition (A ~ lAt for JAF=0-06t). This decrease 
occurs because the spin reduction factor is smaller, and 
the relative electron phonon coupling larger, in the PM 
phase than in the FM phase. For values of A larger 
than the critical value, but not too strong, the mode Q2 
is small and the distorted octahedra only trap a small 
amount of electronic charge, see Fig.®. Interestingly 
the disorder in the position and orientation of the Mn 
octahedra creates states in the energy gap region of the 
insulating ordered phase that favor transport of charge 
making the system to behave as a metal 37, t38j . In 
Figl3a) we plot the density of states for a realization 
of the disordered octahedra phase for X=1.6t in a system 
containing 24x24 Mn ions. The amplitude of the distor- 
tion that minimizes the energy is Q2—0-95. There is not 
gap in the DOS, and there is a finite number of disor- 
der induced states at the Fermi energy. It is interesting 
to compare this DOS with that of the ordered phases 
calculated with the same A, Fig. Elc)-(d). The energy 
gap of the ordered phase disappears when disorder in the 
positions and orientation of the octahedra is introduced. 
The states at the Fermi energy transport electric charge 
and the metallic behavior of this phase is reflected in the 
low frequency spectral weight of the optical conductivity 
cr(w), FiglH^a), that shows a clear Drude-hke character. 

By increasing A the octahedra distortion increases, the 
electronic charge becomes strongly localized, and the sys- 
tem is a disordered mixing of Mn^+ and Mn'*+ . We name 
this state polaronic phase. In FigO a typical octahedra 
distribution and electron charge in the polaronic regime, 
X=2t, is plotted. Note that the electron charge is prac- 
tically localized on Mn ions surrounded by a distorted 
oxygen octahedron. In this strong coupling regime the 
DOS shows a fully occupied polaron band an a strong 
suppression of states at the Fermi energy as can be ob- 
served in Fig. Elb). This is an indication of the insulator 
behaviour of this phase, being more evident in the opti- 
cal conductivity Fig|S{b), that shows a strong reduction 
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FIG. 7: Density of states of the disordered octahedra phase in 
the case of (a) X=1.6t and (b) X=2t. The dased lines indicate 
the position of the Fermi energy. The density of states is 
obtained by solving numerically a system with 24x24 Mn ions 
per unit cell. 



FIG. 8: Optical conductivity of the disordered octahedra 
phase in the case of (a) A=1.6t and (b) X=2t. The optical 
conductivity is obtained by solving numerically a system with 
24x24 Mn ions per unit cell. 



of the low energy weight and the development of a gap 
structure. 

Although is very difficuh to calculate numerically the 
critical value of A for the metal insulating transition in 
the disordered octahedra phase, by comparing qualita- 
tively the DOS and the optical conductivity we estimate, 
A ^1.5t, as the critical value of the electron phonon cou- 
pling (FigPl. 

IV. SUMMARY 

We have calculated the phase diagram of manganites 
at half doping. By combining a realistic microscopic elec- 
tronic model and a mean field treatment of thermal ef- 
fects, we have obtained the free energy of the different 
possible states. With these ingredients we build up phase 
diagrams as function of temperature and electron phonon 
and as function of temperature and antiferromagnetic 
coupling. By changing the value of the electron phonon 
or the antiferromagnetic coupling we simulate different 
manganites with different controlled bandwidth. The 
calculated phase diagrams show a big variety of exotic 
phases some of them, as the CE-CO-00 and the FM- 
CO-00, have been observed experimentally. 



We have obtained that the transition from the high 
temperature paramagnetic phase to the low temperature 
CE phase, depends on the parameters A and Jaf- For 
large values of A the transition can be through a FM-CO- 
00 phase (small Jaf) or trough a PM-CO-00 phase 
(large Jaf)- However for small values of A the transition 
can be direct (large Jaf) or trough a disorder FM phase 
(small Jaf)- 

We have also discussed the nature of the high tem- 
perature PM phase in the regime of finite A. We have 
obtained that in the weak coupling regime the oxygen oc- 
tahedra surrounding the Mn ions are slightly distorted 
and only trap a small amount of electronic charge, be- 
ing the system metallic. However in the strong coupling 
regime half of the octahedra are strongly distorted and 
the charge is localized in the Mn's ions encircled by a 
distorted octahedron, forming a polaron and making the 
system insulator. 
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FIG. 9: (Color online) Schematic view of the disordered octa- 
hedra phase for large values of A. The horizontal and vertical 
lines on the lattice sites represent the positive and negative 
oxygen octahedra distortions. The size of the dots is pro- 
portional to the charge located on the Mn ions. For A=2t 
the octahedra distortion is Q2=H-S5 t, and the electrons are 
practically fully localized in the Mn ions surrounded by dis- 
torted octahedra. 
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